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Abstract 

We derive an expression for the superfluid density of a uniform two-component Fermi gas through 
the BCS-BEC crossover in terms of the thermodynamic potential in the presence of an imposed 
superfluid flow. Treating the pairing fluctuations in a Gaussian approximation following the ap- 
proach of Nozieres and Schmitt-Rink, we use this definition of ps to obtain an explicit result which 
is valid at finite temperatures and over the full BCS-BEC crossover. It is crucial that the BCS 
gap A, the chemical potential /U, and all include the effect of fluctuations at the same level in 
a self-consistent manner. We show that the normal fluid density pn = n — ps naturally separates 
into a sum of contributions from Fermi BCS quasiparticles {Pn ) and Bose collective modes (p^)- 
The expression for p!^ is just Landau's formula for a BCS Fermi superfluid but now calculated over 
the BCS-BEC crossover. The expression for the Bose contribution p^ is more complicated and 
only reduces to Landau's formula for a Bose superfluid in the extreme BEC limit, where all the 
fermions have formed stable Bose pairs and the Bogoliubov excitations of the associated molecu- 
lar Bose condensate are undamped. In a companion paper, we present numerical calculations of 
Ps using an expression equivalent to the one derived in this paper, over the BCS-BEC crossover, 
including unitarity, and at finite temperatures. 

PACS numbers: 03.75.Kk, 03.75.Ss 
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I. INTRODUCTION 

The superfluid density ps is a fundamental signature in all superfluid systems It 
describes the part of the system which does not respond to an external rotation (transverse 
perturbation). Moreover, it is an essential parameter which enters in the two- fluid hydro- 
dynamics of a superfluid, as first discussed by Landau in 1941 j2|. The superfluid density 
ps is quite different from the condensate density n^. In particular, it can be shown that at 
T = 0, the entire system is superfluid {ps = n), in stark contrast to ric 

In this paper, we define and derive an expression for the superfluid density for a two- 
component Fermi gas in the BCS-BEC crossover region at finite temperatures. Our analysis 
is limited to a uniform gas. The calculation of ps is based on the Leggett mean-field BCS 
model of the BCS-BEC crossover, extended to include the effects of pairing fluctuations as- 
sociated with the dynamics of the bound states using the approach of Nozieres and Schmitt- 
Rink (NSR) |4,] . The NSR approximation has also been used to calculate the thermodynamic 
properties in the BCS-BEC crossover at both T = and finite temperatures. As shown in 
detail by Hu et al. jsf 0] , this approximation gives results that are in good agreement with 
Quantum Monte Carlo calculations [3, 0]- Their work gives us confidence in using the NSR 
approximation to calculate the superfluid density at finite temperatures in the BCS-BEC 
crossover, apart from a small region just below T^. where the fluctuations require a more 
careful treatment. 

We note that in the superfluid involved with the BCS-BEC crossover, ps will always refer 
to the number of fermions which participate in the superfluid motion, not the number of 
Bose pairs. Thus at T = 0, = n, where n is the number density of fermions and hence, 
mps is the total mass of the system. 

In the present paper, we define the superfluid density by imposing a "phase twist" on the 
Cooper pair order parameter, endowing the Cooper pair condensate with a finite superfluid 
velocity v^. Following the approach of Ref. ^5 Ps is formulated in terms of the second 
derivative of the thermodynamic potential of the superfluid with respect to Vs- We show 
that the normal fluid density pn = n — p^ naturally separates into a sum of a Fermi quasi- 
particle contribution arising from the standard BCS static mean-field approximation plus a 
Bose contribution arising from the dynamics of the pairing order parameter. The latter con- 
tribution is treated within a Gaussian approximation for the fluctuations around the static 
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BCS order parameter describing the Bose-condensed pair (Cooper) states 10,[ll|. We use 
a single-channel model appropriate for a broad Feshbach resonance Q], which means that 
one deals with an interacting Fermi gas with a varying s-wave scattering length a^. Apart 
from this, our microscopic model is identical to that used in earlier work on the collective 
modes in the BCS-BEC crossover at finite temperatures 

As noted earlier, it is important to keep in mind the distinction between the superfiuid 
density ps and the condensate density ric, related to the average occupancy of the Cooper 
pair state. Numerical results for both and Og are presented in our companion paper as a 

n 

function of both T and [13]. We note that ric has recently been calculated at T = in 
the BCS-BEC crossover using a Quantum Monte Carlo simulation jl^ . 

To date, pg in a Fermi superfiuid with a Feshbach resonance has not been calculated 
by such an ab-initio method, although this has been done for superfiuid ^He as a func- 
tion of temperature 0, [l^. While it is a fundamental property of superfiuids, there have 
been no experimental measurements of ps in the BCS-BEC crossover. As we discuss briefly 
in the concluding section, ps plays a crucial role in two-fluid hydrodynamics. This coUi- 
sional domain should be accessible at finite temperatures near unitarity in the BCS-BEC 
crossover 

We also discuss the equivalence of different formal expressions and definitions for the 
superfiuid density within a given microscopic model. We argue that relating the normal 
fluid density p„ to the thermodynamic potential fi(vs) in the presence of a finite superfiuid 
flow Vs gives a very elegant way of separating out the Fermi BCS quasiparticle contribution 

and the Bose collective mode contribution arising from pairing fluctuations. When 
expanded out in terms of products of single-particle BCS Green's functions (see Appendix 
B), our expression for p^ is extremely complex and not physically transparent. 

We show that the Fermi contribution p^ to the normal fluid density is always given by 
the well-known Landau formula in terms of Fermi BCS quasiparticle excitations. Only the 
values of A and p appearing in the energy spectrum of these excitations change as one sweeps 
through the BCS-BEC crossover. In the BEC limit, the fact that p is large and negative 
means that the Fermi quasiparticles are frozen out by a large effective energy gap over the 
relevant temperature scale ksT ~ /c^Tc <^ |p| and consequently, p^ becomes negligible. 

In contrast, the Bose fluctuation contribution p^ to the normal fluid density becomes 
increasingly dominant as we go from the BCS region to the BEC region, where the dynam- 
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ics of tightly-bound molecules dominate the thermodynamics. Far into the BEC region, 
the Bose fluctuations reduce to the usual Bogoliubov excitations calculated in the Popov 
approximation, which allows for a thermal depletion of the condensate density ^^(T). We 
show in detail how our general expression for the Bose fluctuation contribution to the 
normal fluid density reduces in the extreme BEC limit {\p\ ^ ksTc) to the expected Landau 
expression ^ for the normal fluid density in terms of undamped Bogoliubov-Popov exci- 



tations. This reduction in the BEC limit has recently been proven by Andrenacci, Fieri, 



nation of an expression for defined 
However, we find that there are 



and Strinati .17] based on a direct diagrammatic eva^ 
in terms of a transverse velocity response function 
additional terms in our expression for p„ which are not included in the diagrammatic anal- 
ysis of Ref. llTi]. These terms are negligible in the extreme BEC limit of strongly-bound 
pair states, but become important closer to unitarity where the s-wave scattering length 
becomes very large. 

The present paper concentrates on the formal definition of the superfluid density Ps 
and the derivation of an explicit (but still formal) expression for a specific microscopic 
model which includes contributions from the Fermi BCS quasiparticles and the Bose pairing 
fluctuations. We concentrate on the structure of these two contributions to p„ and the 
underlying physics of the pairing fluctuations which give rise to p^. A companion paper by 
the authors presents the results of extensive numerical calculations of our expression 
for function of both the temperature and s-wave scattering length a^. In such 

calculations, it is important to use the renormalized values of BCS gap A and the Fermi 
chemical potential p within a Gaussian approximation that includes the effects of the same 
pairing fluctuations which describe the Bose collective mode contribution p^ to the normal 
fluid density. 

II. FORMAL EXPRESSION FOR THE SUPERFLUID DENSITY 

Our expression for the superfluid density is based on the equilibrium thermodynamic 
potential for a current-carrying superfluid. Thus, our starting point is the partition function 

Z = f V[i:,ij]e-^^^'^^ (1) 
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expressed as a functional integral over fermionic Grassmann fields if) and if) jiol . 
imaginary-time action 5'[?/',?/'] is given by 



dr 



j dr ^'ilj^{x)dripa{x) + H 



The 



(2) 



where /? = l/ksT. Here, we use the notation x = (r,r) where r denotes spatial coordinates 
and T = it is the imaginary time variable. We set h = 1 throughout. H is the usual BCS 
pairing Hamiltonian, 



H = J dr Y.^<^(^) (2!^ ^'"'^^^ ~ ^ / 



(3) 



U is the parameter characterizing the s-wave scattering interaction between fermions in 
the two different hyperfine states, denoted by the spin indices a =t, i- From the Lippman- 
Schwinger equation for the two-body scattering problem, U is related to the s-wave scattering 
length tts by 



1 

U 



m 
Ana. 



-1 



(4) 



where = k^/2m. Throughout this paper, we take the volume to be unity. Our analysis 
is restricted to uniform gases. 

The Bose pairing field A{x) that includes fluctuations about the mean-field static BCS 
order parameter A is introduced through the Hubbard- Stratonovich transformation. 



'l?[A, A*]exp 



/3 



dr / dr 



\A{x) 



A*{x)iJi^^ + A{x)ij^ij^ 



(5) 



With this identity, the partition function becomes 



J V[ij, ij]V[A, A*] exp j^T j dr ^ ^.(x) {d^ + ^- j ^,{x) 



\A{x) 



U 



j V[i),i)\D[A,A*]e-KY> 



|A(x) 
U 



(6) 
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where we have introduced the Nambu spinors 

/ 



and is the inverse of the 2x2 matrix Nambu-Gorkov BCS Green's function, 



(7) 



V 



G-Hx,x') 



6{x-x'). (8) 



The integration over the Grassmann fields ip in Eq. ® can be performed in straightfor- 
ward fashion to give 

Z = J P[A, A*]e-^««, (9) 



where ^ 

^efr[A, A*] = dr J dr - Trln[-G-i]. (10) 

The trace in Eq. (fTn|l is performed over space and imaginary time variables, in addition to 
the Nambu indices. We have used the standard identity InDetA = Trln A. 

The key function of interest in this paper is the thermodynamic potential fl, defined by 

n = -kBT\nZ. (11) 

All thermodynamic quantities of interest can be calculated once Q is given in some mi- 
croscopic approximation. We make use of the fact that ps can also be obtained from the 
thermodynamic potential of a current-carrying superfluid. To impose a current, one applies 
a "phase twist" 9] to the order parameter A(a;): 

A(x) A(x)e*^-^ (12) 

The superfluid velocity associated with this imposed phase twist is 

= — , (13) 

where M = 2m is the Cooper-pair mass. Treating Q as small, the superfluid density is 
obtained from the lowest-order change in the free energy of the system {F = Q + fiN) due 
to the added kinetic energy of the imposed superfluid flow j^. This extra kinetic energy is 



Af^f(Q)-F(0).f (^) (14) 



with 



Note that the superfluid density defined here is the superfiuid number density and not the 
superfiuid mass density used in discussions of two-fiuid hydrodynamics. As can be seen from 
Eq. ()14|). psin is the total mass involved in the superfiuid fiow, with m being the Fermi atom 
mass. 

For a fixed number of fermions A^, 

Microscopically, Q can be expressed as a functional of the mean-field gap A, the chemical 
potential /i, and the phase twist Q. In addition to an explicit Q- dependence, Q also depends 
on the phase twist implicitly through the gap A{Q) and the chemical potential /u((5). Using 
these facts, we can write Eq. (fT^ as 



dQ^Q^, V^^?Va,^ \9^J,9Q' Kdy^J^dQ^ dQ' 



'd^A 








m 





(17) 



In going from the first to the second line of Eq. (jl7|) . we have made use of the number 
equation 

fdn\ fdn\ fdn\ fdA\ 

The evaluation of the derivatives at Q = is left implicit on the right-hand side of Eq. ()17p 
and we have also made use of the fact that the first-order corrections to /i and A vanish: 
{dfi/dQ)Q^o = {dA/dQ)Q^o = 0. Separating the mean-field and fiuctuation contributions, 
Q = Q^{+6Q, where = Q{A{x) — » A), we can write the derivative of the thermodynamic 
potential with respect to A as 

The first term in Eq. (jl9|) vanishes since, by definition, dQ{A{x) A)/dA = 0. As argued 
in Ref. the second term in Eq. ()19|) is a higher-order correction, beyond the Gaussian 
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theory we use to evaluate fl in Section III. Thus, for the sake of consistency we ignore this 
contribution and set the second term in Eq. (|T7|) equal to zero. Our final expression for the 
superfluid density is thus 

In Eq. ()20j) and elsewhere, the constancy of A and fi in taking derivatives with respect to 
Q is left implicit. This formula is the basis for our discussion of ps in this paper. 

We note that by ignoring terms proportional to (dQ/dA) in Eqs. ()17|) and ()18|). the 
number equation we use to define the superfluid density reduces to 

This expression keeps A fixed, meaning that derivatives of the form (dA/dfi) do not en- 
ter into the resulting equation for A^, in contrast to Eq. (jl8p . which includes higher order 
corrections. For our calculations to be consistent, the chemical potential used to evaluate 
our expression for ps must also be calculated using Eq. as done in Refs. [lOllllI- The 
contribution of the higher order term {dQ/dA) to the calculation of equilibrium thermody- 
nan,ic quantMe. ha. been di.cu.ed »n,e .ecento^e. gflQ. In particular in the 
context of the BCS-BEC crossover problem, Refs. p, m{ make use of the full number equa- 
tion given by Eq. (fTHj) to obtain results that are in excellent agreement with Monte Carlo 
simulations at both T = [7J and finite T We defer further remarks on this to Section 
VI. However, it appears from Ref. j6| that this derivative brings in the effect of cubic and 
quartic fluctuations 221 which have the effect of renormalizing the strength of the effective 
interaction between stable Cooper pairs |23| . 

In Appendix A, we review the arguments demonstrating the equivalence of Eg. ([2()ll and 
the usual definition of Ps in terms of the transverse current correlation function 18|. 



III. THERMODYNAMIC POTENTIAL FOR A CURRENT-CARRYING SUPER- 
FLUID 

In order to calculate the thermodynamic potential for a current-carrying superfluid, the 
phase twist is applied to the order parameter that enters the inverse Green's function 
in the action given by Eq. (fTUI). To remove the phase from the order parameter, we apply 
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the unitary transformation G ^ = U 2^ 3, H]' where 



I- (22) 

Owing to the invariance of Tr ln[— G^^] with respect to the action of a unitary transformation 
of G~^, the effective action with a phase-twisted order parameter can be written as 

^eff[A,A*,Q] = I (ir ^^-Trln[-G-^], (23) 

where (p = — iV) 



G-\x,x')=\ 2m ^ U(x-x'). 24 

The thermodynamic potential for a current-earring superfluid can be evaluated from this 
action, using Eqs. Q and (fTTj) . once some approximation is introduced so that the functional 
integration in Eq. © can be carried out. Following the standard prescription, we expand 
the action in powers of fluctuation about the mean-field BCS pairing field: A(x) = A + A(x); 
G-i = Go^ + S, where Gq^ = G'\A{x) A) and 

/ 



A{x) 
A(x) 



6{x-x'). (25) 



Clearly A(a;) corresponds to the fermionic self-energies due to coupling to Bose collective 
modes involving pair fluctuations in the Cooper pair channel. 

Using the expansion Trln[-G~i] = Trln[-Go^(l + GqS)] = Trln[-Go^] + Trln[l + 
GoS] = Trln[-Go^] + En=i Tr[(GoS)"](-l)"+V^. we expand Eq. ^ up to quadratic 
order in the Bose fluctuation field A to obtain the Gaussian action, S'causs = S^'^^ + S^'^\ 
Fourier-transforming, the mean- field S^^^ and fluctuation S^"^^ contributions are given by 

^(°) = /5^-5:trln[-Go^(fc)] (26) 
^ k 



and 



^(2) ^ + Vtr[Go(A:)S(-g)Go(A: + g)S(g)] 



k ^ 2 /j q 

VAtMA. (27) 

q 



In Eq. ()27p . q = (q, zz/m) and k = (k, icj^) are 4- vectors denoting the momenta q and k 
as well as the Bose and Fermi Matsubara frequencies ium = 2Tim/ (3 and iuJn = 27r(n + 
1)//?, respectively, where m,n are integers. In momentum-frequency space, the mean-field 
(denoted by the subscript "0") Nambu-Gorkov BCS Green's function Go(^) foi' the current- 
carrying BCS superfluid is defined by its inverse, 

G-W.(,:.„-!i^)-(a + |l)f3 + Af. (28) 

Here, = k^/2m — fi, while Ti,f3 are Pauli spin matrices. We have assumed that the 
mean-field order parameter A = A* is real. In the last line of Eq. ()27p. we have defined 
the spinor = (A(g), A(— g)), and the matrix elements of the inverse 2x2 matrix pair 
fluctuation propagator M for a current-carrying superfluid are given by 

Mii(g) M22(-g) 11^^ // , A/^ fi\ fnn\ 
^ = ^ = ^ + ■^2^G'o,ii(/i; + g)Go,22(A;) (29) 

and 

Muiq) M2i(g) ^^r' (i , \n (i\ 

- ^ -3 2^Go,u{k + q)Go,i2{k). (30) 



(3 P P 



k 



Here, Go,jj denotes the ij-th element of the matrix mean-field BCS Green's function defined 
by Eq. 

Substituting Scs ~ S^'^^ + S**^^^ into Eq. Q and performing the Gaussian integration 
over the Bose fluctuation fields (A, A), the thermodynamic potential for a current-carrying 
superfluid reduces to 

n{Q) = ^-^EtrM-Go^(A:)] + ^ElndetM(g) 

= fi^(Q) +1]^(Q). (31) 

This formula will be used to calculate ps in Eq. (j2(Jj) and thus plays a key role in the rest of 
this paper. The first two terms of Eq. (|!TT|l comprise the mean-field contribution from Fermi 
BCS quasiparticles, 

Q^iQ) = ^-lT.^rln[-G,\k)]. (32) 

1"^ k 



We emphasize that the values of A and fx in Eq. ()28|) evaluated using our Gaussian theory are 
strongly renormalized from their mean-field values by the effects of fiuctuations in the Cooper 
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pair field and the associated collective modes, as described by the NSR theory [lOl, 
The values of these microscopic parameters for a current-carrying superfluid are obtained 
by self-consistently solving the gap equation, (dS^^ydA) = 0, together with the number 
equation N = -(l^(Q)/9/i)A, where S^^'^ is given by Eq. and i7(Q) is given by Eq. (j31|) . 
Recall that our expression for the superfluid density, given by Eq. ((201), leaves A and fj, 
fixed, so we only require the values of these quantities in the current-free state, found from 
a^WfO = 0)/gA = and = -{n{0)/dfi)A- Further details of this calculation are given 
in Ref. [13|. 

The contribution from the Bose collective modes in Eq. PT| is 

^^''(Q) = ^ E lndetM(q,2i/™). (33) 

The collective modes will be shown to play an increasingly important role in p„ as one goes 
from the BCS to the EEC regimes. The spectrum ujq of the collective modes is determined 
from 

detM(q, tu^ ^00^ + «0+) = Mu{q)Mu{-q) - MUq) = 0, (34) 

where — > cUq + iO~^ denotes the usual analytic continuation from imaginary Bose frequen- 
cies. In most of the BCS-BEC crossover, these collective modes will be damped at finite 
temperatures (i.e., cjq has an imaginary part). 

IV. SUPERFLUID DENSITY IN THE BCS-BEC CROSSOVER 

In this section, we derive an explicit expression for the superfluid density in the crossover 
starting from the formula given by Eq. ()20j) for the model defined in Section III. From 
Eqs. fl28|) and (jHTj) . one sees that the thermodynamic potential for a superfluid with a finite 
superfluid velocity = Q/M is equivalent to the thermodynamic potential for a current- 
free superfluid (v^ = 0), but where the chemical potential and Matsubara frequencies are 
now Doppler-shifted j25|: 

/i — s> /i — Q'^ /8m = jl, (35) 
iujn i^n — k ■ Q/2m = 2u)„, (36) 
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iPm ^ium-q - Q/2m = iu.m. (37) 

Considering separately the effects of the shifts to the chemical potential and the Matsubara 
frequencies, we can write the second-order derivative of Q with respect to Q (keeping A and 
yU fixed) as 

1 dn Q d'n ^ ^gg^ 



4m (9/i 2m djldQ \ dQ"^ 



Evaluated at Q = 0, the middle term in Eq. (jHHj) vanishes and Eq. reduces to 

In the last line, we have made use of the number equation n = — (<9f2/<9/i)^,A.Q-»o = 
— {dQ/dfi)A- Since n = ps + pn, Eq. gives us the following expression for the nor- 
mal fluid density: 

""^ 

Carrying out the summation over Fermi Matsubara frequencies in Eq. fl32|) . the mean- 
field BCS quasiparticle contribution to the thermodynamic potential in the presence of a 
current becomes j^l 

fi^(Q) = ^ + E (a - ^k) - I E In [l + e-^('^-Q/2™+^'')] , (41) 



where the single-particle quasiparticle energies are given by = \/^k + ^^ with = 
k^/2m — p, where p is the Doppler-shifted chemical potential defined in Eq. ()35|). 

Summing over the fermion Matsubara frequencies in Eqs. ()29|1 and ()3()|1 . the matrix ele- 
ments of the inverse matrix propagator for pair fluctuations in the current-carrying superfiuid 
are given by 



Mii(g) = M22(-g) = ^ + E 

^ k 



fk f] 



+ /k - /] 



k+q 



k+q 



'"k'^k+q 



^'k^k+q 



ium-q - Q/2m + -Ek + -Ek+q 



Wm-fl - Q/2m - -Ek - -Ek+q 
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/k+q) 



/k+q) 



'"k'^k+q 



q • Q/2m + -Ek - -Ek+q 



"k^k+q 



iz/m - q ■ Q/2m - Ek + -Ek+c 



(42) 



and 



Mi2(g) = M2i(g) = E 



k+q 



Wk^^k^^k+q^k+q 



+ (/k+q /k ) . 



ij^m - q ■ Q/2m + E'k + -£^k+q 

^^k^^k^^k+q^^k+q 



+ (/k^ /k+q) 
+ (/k ~ /k+q 



zz/^ - q ■ Q/2m - -Ek - -Ek+q 

^k^k^k+q^k+q 

zz/„ - q ■ Q/2m + ^k - ^k+q 

^k^k^k+q^k+q 

ium - q ■ Q/2m - Ek + -Ek+q. 



where 



/p±^/(p-Q/2m±E, 



(43) 



(44) 



are the Fermi distribution functions. Here, Up = y {1 + ^p/Ep)/2 and f p = y (1 — ^p/Ep)/2 
are the usual Bogohubov quasiparticle amphtudes. Recall that the normal fluid density is 
evaluated at fixed fi and consequently, the dependence of on Q can be ignored for the 
sake of calculating p„ in Eq. (jlU)). The expressions given by Eqs. and pHj) reduce to 
the standard expressions in the literature |lO| for Mjj(q, iz/m) when = 0. 

The distribution functions appearing in Eqs. ()42|) and (pSj) involve Doppler-shifted Fermi 
quasiparticle energies: p ■ Q/2m ± Ep. The shift p ■ Q/2m reflects the fact that additional 
Fermi quasiparticles will be excited when the superfiuid velocity is finite since thermal 
equilibrium is defined with respect to the stationary lab frame j25|. 

Using the thermodynamic potential in Eq. (j3H) . the normal fluid density pn is given by 
the sum of Fermi quasiparticle and Bose collective mode contributions: 



where 



Pn 



m ^ 



Pn+Pn, 



k Q 



m 



tr[Go(A;)Go(fc)] 



and 



Pn 



2m ^ 1 
q (detM) 



~ \ 2 



detM 



detM\ 



'(9detMy 



(45) 



(46) 



(47) 
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Here, Q = Q/|Q|. The expression for the Bose contribution is very compactly given in 
terms of the determinant of the inverse fluctuation propagator detM, the zeros of which give 
the spectrum of the Bose collective modes. The simplicity of this expression for is lost 
when expanded in terms of products of current-free BCS Green's function (see Appendix 
B). One can show after a little work that the result given by Eq. (|47p is identical to that 
obtained in Ref. based on a calculation of the current response to a superfiuid flow. 

The normal fluid density due to Fermi BCS quasiparticles given in Eq. ()46|) is readily 
identified as the long-wavelength, static limit {q — > 0) of the BCS current-current correlation 
function (multiplied by —m). Carrying out the Matsubara frequency sum in the usual way, 
Eq. (pUj) reduces to 



Pn = 77rTT(k-Q) 



mJ (27r)3' dEk 
±f^y(yjm], (48) 



3m 7 (27r)3 \^ dEk 
This is the well-known Landau formula for the normal fluid density of a uniform weak- 
coupling BCS superfiuid, arising from thermally-excited Fermi BCS quasiparticles In 
our case, it is valid for the entire BCS-BEC crossover, taking into account that the quasi- 
particle spectrum depends on A and p which are renormalized from their mean-field BCS 
values by the inclusion of the effects of Bose fluctuations 000 . Note that the Landau 
formula given by Eq. pHjl also results by using Eq. (jHj) in Eq. (Hn|) . 

Eq. (jTTj) describes the contributions to the normal fluid density from fluctuations 6 A of the 
Bose pairing field. In general, the Bose pair excitations are damped at finite temperatures, 
coupling to the continuum of BCS quasiparticle states. As a result, the Bose fluctuations 
will have a finite lifetime and p^ will not reduce to the usual Landau formula involving Bose 
excitations. In the BEC limit, however, the pair binding energy becomes very large and 
BCS quasiparticles are strongly suppressed. As a result, damping will not occur. In this 
limit, we expect that our expression for p^ will be given by Landau's formula for a Bose 
superfiuid. In the next section, we give the details of this proof. 



V. THE NORMAL FLUID DENSITY IN THE BEC LIMIT 



Close to unitarity and on the BCS side of the BCS-BEC crossover. Landau damping of the 
Bose collective modes described in Section IV arises due to scattering processes that involve 
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BCS Fermi quasiparticles: cDq+E'k = -Ek+q. When such damping occurs, the collective modes 
strongly hybridize with BCS Fermi quasiparticles and the concept of well-defined, long-lived 
Bose excitations breaks down. In this region, we do not expect the normal fluid density 
to be given by a Landau formula for the Bose excitations. In the strong-coupling BEC limit. 



however, the fermions form bound pairs with a large binding energy j20|: -^binding = — 
As a result of this large binding energy, Fermi quasiparticle excitations, which involve the 
breakup of pairs, become completely frozen out over the experimentally relevant temperature 
scale ksT ~ ksTc <^ -^binding- Because the Fermi quasiparticles are frozen out, they no longer 
contribute to Landau damping of the Bose collective modes. Thus, in the BEC limit, the 
normal fluid is comprised of a gas of well-deflned Bose excitations and one expects p,f will 
reduce to the usual Landau expression for Bose excitations in this limit 2^. In this section, 
we show how this result emerges from our formalism (which is valid in the entire BCS-BEC 
crossover) in the BEC limit. 

Deep in the BEC region, the chemical potential becomes increasingly large and negative. 
In the strong-coupling limit where A, Q^/Sm, /c^T ^ the BCS gap equation for the 
current- carrying superfluid {Q = Mvg), 

^ = lT.GoMk), (49) 

i"^ k 

can be solved analytically. This gives jl = p — Q'^/8m = — l/(2ma^), which is one-half the 
molecular binding energy Q]. When ksT, the BCS quasiparticles are frozen out 

{fp — ^ 0; /p ^ 1) and the BCS quasiparticle contribution p^, given by Eq. (jlHj) . vanishes. 

In the low-energy regime <^ the spectrum of Bose excitations is expected to have 
the form y^(cq)^ + (q^/2m*)^. To extract the contribution of these modes to the normal 
fluid density p^, we set /+ = and = 1 in Eqs. (jSH)) and (jHO)) and then expand the 
inverse fluctuation propagator matrix elements in powers of q. This procedure gives 

MllM. + 5|q|2 + C{iiy^ - q ■ Q/2mf + D{iUm - q ■ Q/2m), (50) 

and 

~ A + F|q|2 + G{iu^ - q ■ Q/2m)^. (51) 

We note that outside the BEC region, where ksT ~ C(|/i|), we cannot set /+ = 0, = 
1. Consequently, the terms in the inverse fluctuation propagator responsible for Landau 
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damping, given by the last two lines in Eq. ()42p and Eq. ()43|) . cannot be neglected. In this 
case, it is well-known that one cannot carry out an expansion in powers of q and iu^, as 
in Eqs. (I5UI) and (|51|). since these terms are singular in the long wavelength, zero frequency 
limit [30, [sil. This means that the expansions in Eqs. (j5(J|) and (j51|) are not valid in the 
unitarity or BCS regions. 

Apart from the shift to the chemical potential given by Eq. (jHSj) . the expansion coefficients 
in Eqs. fl50|) and ()51|) are the same as for the Q = case given in Ref. 10|, namely 

A2 



(52) 



5 = E 



k L 



2-3^ I — + 



El 



m 



kP cos^ ' 



A2 

-2 + 13^-10^ 
El 



k/ J 



16^3' 



(53) 



(54) 



D 



E 



k 4^3' 



(55) 



E 



El ^ 



IkP cos^ 1 









V El 


Ei \ 



IQEl' 



(56) 



and 



G 



IQEl 



(57) 



We have made use of the gap equation, given by Eq. to eliminate l/U from (/5)-^Mii(g). 
In the strong-coupling BEC limit, A ^ and we can further expand the integrands in 
powers of A/|/i|. To leading order, using = (2ma^)~"'^, we find 



k 



IGtt 



(58) 



B 



E 

k 



8m^l 



|kp COS^ : 

4m2|k 



mas 
32^ 



(59) 



-E^ 

k 



1 m^a? 



IGtt 



(60) 
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and 



n V- 1 ^^^^ 

D^-i:—,=-^- (61) 



k 



To leading order, we find F ~ A^a^ and G ~ A^aJ, which are vanishingly small in the BEG 
limit, as — * 0. Similarly, since C oc (ma^)^, we set this coefficient equal to zero as well. 
However, since A^ oc a~^, one finds that A oc a^, and we retain A in Eqs. (I5(J|) and (j51|) . 
With coefficients given by Eqs. (I^Uj). and (|HT|). and setting C = F = G = 0, we find 

M 



detM(q, iV^) = 2ABq2 + E^q^ - (tu^ - ^ 1 . (62) 



Since the fiuctuations spectrum is given by the zeros of detM(q, tUq), one finds 



l2AB „ 52 



Wq(Vs) = q ■ V, + y -^q2 + — q4 



q- + 



We note that the value of A appearing in this expression is temperature-dependent. Since 
A(T) 7^ is associated with the existence of a molecular Bose condensate in the BEG region 
of interest, the dispersion of Bose collective modes can be written in terms of the condensate 
density ric. In Ref. |13|], we show that the corrections 6nc to the mean-field expression for 



the condensate density. 



^co{T) = E tanh2(/3Ek/2), (64) 



are negligible throughout the BGS-BEG crossover, within our NSR Gaussian approximation. 
Thus, we can use Eq. ()64|1 to determine the condensate density in the BEG limit (where 

n.iT) = . (65) 

It is important to emphasize that in obtaining this expression, we have only taken the 
limit \^\/kBT 00, where tanh.^((3Ek/2) 1. However, A still has a strong temperature 
dependence arising from the thermally excited pairing fiuctuations which are not frozen 
out. The temperature dependence of A(T) is calculated within our Gaussian approximation 
throughout the BGS-BEG crossover in Ref. jisj ]. 
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Using the result in Eq. ()65|) . one can show that the sound velocity in Eq. ()63|) can be 
written as 

^2 ^ A(T)a^ ^ UMn.jT) 

AM ^ ^ 

for an interacting gas of bosons of mass M = 2m. This is the standard Bogoliubov-Popov 
sound velocity with Um = AnaM/M, but with the molecular scattering length given by the 



mean-field result a^/ = 2as 



20|. 



In order to get the correct value of the molecular scattering length qm — O.Gag in the EEC 
limit 2^, one would have to include the effects of 4-body correlations which are beyond the 
2-body physics contained in our Gaussian theory; i.e., we would need to expand the action 
to quartic order in fluctuations As pointed out by Hu et al. ^, the correct renormalized 
value of ajvf emerges when one calculates /i using the number equation given in Eq. (fTHj) that 
includes the contribution from dfl/dA. Thus, while we do not consider it in this paper, it 
appears that we understand how our present calculation can be improved to get the correct 
value of aM — CGa^. 

Using the expression for detM(q, zi>„J given by Eq. (j62p . it is straightforward to evaluate 
in Eq. (gZI). Making use of Eq. ^,Eq. ^ reduces to 



detM(q,ii^m) 
Using this expression, we find 



{iUrn - q • Q/M)^ - u^ivs 



0) 



(67) 



1 



detMl ^ 



detM 







2 










q Q 

M 




(68) 



where = uj^{Q = 0) is the usual Bogoliubov-Popov excitation energy in the absence of a 
superfiuid flow, = 0. Using this result in Eq. (jTTj) . and recalling that vanishes in the 
BEG limit, we obtain 



Pn = Pn 



M / q-Q ^ 

3 ^ \ M 



2iip^Y + 2ujt 



[IVrr 



OJr, 



|2- 



(69) 



To bring out the physics of Eq. ()69|). it can also be written in terms of the transverse current 
correlation function for a dilute Bose gas of interacting molecules |29|. 



Pn 



-J ("3^1 tr [D(q,iz/„)D(q,iz/^)] 



(70) 



18 



where 



D(q,2i/„ 



y ijym- 



(71) 



is the 2x2 Bose propagator describing the Bogohubov excitations. Carrying out the Bose 



frequency sum in Eq. (fTUj) as in Ref. 



29|, we find the expected result 



2 r d^q 2 / dfT'sij^q) 



3MJ (27r)3'l [ duj J' ^^^^ 
where nsiuj) = {e^^ — is the Bose distribution function. Equation ()72|) is precisely 
Landau's formula for the normal fiuid density of a Bose gas described in terms of Bogoliubov 
excitations [2^. Recall from Section I that ps and hence Oj, always refers to the number 
of fermions. Thus, Eq. (f7^ is twice the usual expression [2^, reflecting the fact that it is 
counting the number of fermions (not the number of bosons) involved with a normal fluid 
composed of Bogoliubov excitations of a molecular BEC. 

As we have already discussed, retaining terms in Eq. (fTHj) that are proportional to 
{dVL/d/S) leads to the renormalization of the molecular scattering length, from = '^cls to 
qm — O.Gag. To be consistent, one must include the analogous terms in Eq. ffTTj) and addi- 
tional terms will be generated in our deflnition of the superfluid density given by Eq. ipiH) : 
Ps Ps + 4:m{dQ/dA)f^[d'^A/dQ'^ — {dA/dfi)d'^fi/dQ'^]. In the extreme BEC limit, however, 
d'^A/dQ'^ ^ as the BCS quasiparticles become frozen out, and d'^fi/dQ^ l/Am [as 
shown below Eq. so that ps — > ps + npj^A, where Upj^A = {dQ/dA)^{dA/dp) is the 

correction to the number equation jQ]. Using this new expression in Eq. (jH^ . Upf^^ just 
adds another contribution to the total density n and Eq. (jiUI) remains unchanged in the 
BEC limit. Thus, even if we retain terms in Eqs. (|T7jl and p8|l that lead to the renormal- 
ization of the molecular scattering length, in the BEC limit, p„ is still given by Eq. (001) ■ 
Consequently, our major result in Eq. ()72|) still holds in the BEC limit when we include the 
higher-order corrections, except that will now be ~ O.Qag. 

For completeness, we write down the pair fluctuation contribution to the thermodynamic 
potential in Eq. (j33p in the BEC limit. Using Eq. (j67p . the Bose Matsubara frequency sum 
can be evaluated analytically [2^ and we flnd 

fi^(Q) = ^ E + ^ E In [1 - e-^('^-Q/^^+-<^)] . (73) 

q P q 
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Inserting Eq. (f75jl into Eq. pUj) also leads to Eq. (f?^ . The integrands in Eq. (fTHjl are strictly 
only valid at small momenta, such that q^/2M <^ A, and the divergent zero-point energy 
contribution to Eq. (f7H|) must be regularized by an appropriate choice of cutoff. 

The result given in Eq. (j72j) for the BEC limit of the BCS-BEC expression has also been 
derived in Ref. using a diagrammatic approach. This is discussed in Appendix B. 



VI. CONCLUSIONS 



In the present paper, we have derived an explicit formula for the normal fluid density p„ 
in terms of two contributions. One is the expected contribution given by Eq. ()46|) arising 
from Fermi BCS single-particle excitations. In the BEC limit, vanishes since the effective 
quasiparticle energy gap becomes very large. Physically, the pair binding energy becomes 
very large and the Fermi quasiparticles, which are excitations corresponding to the breakup 
of these pair states, become frozen out. 

The most interesting contribution to p„ in the BCS-BEC crossover is the contribution p^ 
from collective modes associated with the dynamics of the pair states. This is given in our 
NSR formalism by Eq. (^Tj). Within this Gaussian approximation to the pair fluctuation 
propagator, as summarized in Eqs. and (jlHj) . one can proceed to calculate numer- 
ically. The results of such calculations are discussed in a companion paper [l3| over the 
whole BCS-BEC crossover and as a function of temperature. 

The delicate nature of the Cooper pair molecule in the unitarity region of the crossover 
leads to damping of the collective modes given by Eq. (jS^)- This means that, in general, 
p^ is not given by a simple Landau expression such as Eq. (f?^ . However, one does expect 
such a Landau formula to emerge in the extreme BEC limit where the Cooper pairs become 
very strongly bound and the system reduces to a weakly interacting Bose gas of stable 
molecules. In Section V, we showed how this expected result emerges naturally from our 
general formalism. 

In deriving our key starting formula for the superfluid density given by Eq. ()20|) . we 
neglected the contribution of (dfl/dA) in Eqs. (fTTj) and (fTBj) . We argued that such a term is 
a higher order correction which cannot be consistently included in a Gaussian approximation 
on which our formal analysis and numerical calculation are based. The role of the second 
term in Eq. (fTH|) has been discussed in recent calculations 0,0,1211. In particular, Hu, Lui, 
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and Drummond [6| have shown that in the BEC hmit, this term in the number equation 
gives rise to a renormahzation of the molecular scattering length from the mean-field value 
qm = 2cuto the correct value — O.Ga^. This result is consistent with the calculation of 
Ohashi |2j] who went past our NSR Gaussian pairing fluctuation approximation to include 



33, 



the effects of cubic and quartic fluctuations (for a diagrammatic analysis, see Refs. 
Is^). Ohashi found that these higher order effects lead to a renormahzation of the effective 
interaction between molecules, and obtained a value for close to the result of Petrov et 
al. 0- 

We conclude that the neglected contribution to Eq. (fTH|) picks up an important class of 
fluctuations left out of our Gaussian model, which are precisely those needed to give the 
correct molecular scattering length in the BEC limit. As one knows from other problems, 
derivatives of the Gaussian thermodynamic potential can generate results which describe an 
improved model. This emphasizes the usefulness of calculating starting from the result 
in Eq. (j20|l . 

The superfluid density was first introduced by Landau in connection with a two-fluid 
theory for the coUisional hydrodynamics of a Bose superfluid The form of the Landau 
two-fluid hydrodynamics is generic for any superfluid with a two-component order parameter 
(amplitude and phase) yji]. The frequencies of the resulting hydrodynamic modes are given 
completely in terms of the equilibrium thermodynamic functions, including the superfluid 
density. The precise values of these equilibrium quantities depend on the nature of the 
dominant thermal excitations, which can be different in different superfluids. In the BCS- 
BEC crossover, one goes from the BCS limit, where Fermi BCS quasiparticles dominate 
the thermodynamics, to the BEC limit where the Bose collective modes (Bogoliubov-Popov 
excitations) dominate the thermodynamics. 

As a result, a careful discussion of the two-fluid collective modes requires a careful anal- 
ysis of the changing weights of the Fermi and Bose excitations as we pass through the 
BCS-BEC crossover, both for thermodynamic quantities such as the entropy and compress- 
ibility as well as the equilibrium superfluid density. The advantage of calculating ps from 
the second derivative of the thermodynamic potential f2(vs) calculated within a Gaussian 
approximation for the fluctuations, as we do in our work (see also Ref. ^^), is that all other 
thermodynamic functions can also be determined from VL{vs = 0). Heiselberg 3^ has given 
an informative first study of first and second sound in the BCS-BEC crossover for a uni- 
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form gas by calculating ps and other thermodynamic functions in the BEC and BCS limits 
and interpolating into the unitarity region (|as| — > oo). We hope to give a more definitive 
discussion of first and second sound using the numerical results for ps given in Ref. Q. 
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APPENDIX A: SUPERFLUID DENSITY AND THE TRANSVERSE CURRENT 
CORRELATION FUNCTION 



In this section we review the relationship between our definition of the superfiuid density 
and the transverse response definition commonly evoked in the literature Ql • We start with 
the partition function expressed in terms of both Bose and Fermi fields, given by Eq. ®. 
Applying a phase twist to the order parameter as was done in Section III, Eq. ^ is written 
as 



Z = j P[V;,V^]P[A*, A]e-^[^''^'^*'^'Q] 



where 



5[^,^,A*,A,Q] = / d^x 



<i!\x) [-G"^(x,x')l *(x') + 



|A^ 
IT 



and G ^(x, x') is given by Eq. (j2H). From Eq. (j^ . keeping A and p fixed, 

— ' — Ta d[x — x) 



and 



dQ 



d^G-Hx,x') 



2m Am 



Using Eqs. (lA2l) . (jA3jl . and (lA4l) . we obtain the relations 

d'^x ^Hx 



dS_ 

dQ 



2m Am 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 
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and 



(fx 



X] 



1 

4m 



fsj 5{x- x') 



We use these expressions to obtain the relation 



'dSV d^S 



-,-s 



Q=0 



OQ J dQ^ 

— [ d^x *^f3* 
Am J 

J d'^xd'^x' l^\x)%^<if 

( /? 



'x) 



2m 



where 



4m^~i /rf'a;dVj,(a;)j,(x'))e' 



N = I dv $:^,(r)^,(r) 



is the number operator, and, having arbitrarily chosen Q = z. 



is the z-component of the current density operator. 

Using Eq. ()A7|) along with the thermodynamic potential in the presence of a 
flow. 



we obtain 



n{Q) = -Tin / V[tlj,ij]V[A*,A]e-^^^''^'^*^^'^\ 



T 
Z 



-T 



Z 



V[^,mA\A] l^-^N -^jd'x d^x' h{x)h{x') ) e 



23 



where = {N)q and = /3 / d'^x j^. {■ ■ ■)o denotes an expectation value with respect 
to the current-free state: 

(...) = ! /r'[v5,^]p[A*,A](---)e-^[^'^'^*'^''^=°l, (A12) 
ZoJ 

and 

Zo = J V[^, m^*, A]e-^['^'^'^*'^'Q=°1. (A13) 

Note that the second hne in Eq. ()A11|) vanishes by symmetry, i.e., (72)0 = 0. 
Comparing Eq. (jAllj) with Eq. pO|) . we obtain the result 

Ps = n- m{jJ^)o, (A14) 

which identifies jl^ 

Pn = m{J^J^)o (A15) 

as the normal fiuid density. 



APPENDIX B: NORMAL FLUID CONTRIBUTION FROM THE BOSE FLUC- 
TUATIONS 



A central result of our paper is the formal expression in Eq. ()47|) for the Bose fiuctuation 
contribution to the normal fiuid density, where the M matrix elements are given in Eqs. (j42|l 
and In this appendix, we "unpack" this formal result for p^ to give it more explicitly in 
terms of single-particle Nambu-Gorkov Green's functions of a current-free BCS superfiuid. 
This will allow us to compare with other results in the literature [iT^ . 

Expanding the Q-derivatives in Eq. (jlTj) , the Bose fiuctuation contribution to pn is given 

by 



2m ^ 1 
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Here, Mij = Mjj(Q = 0). To express Eq. ()B1|) in terms of current-free Green's functions, 
we use the identities 



dQ 



2m 



[Go,n{p)Go,n{p) + ^0,12(^)^0,12^)] , (B2) 



dQ 



2m 



[GoMp)Go,uiP) + GoMp)GoMp)] > (B3) 



dGoMp) 



dQ 



P Q 

2m 



[Go,22(p)Go,22(p) + Go,12(p)Go,12(p)] , (B4) 



o'GoMp) 



2m 



ii,Q-*o 
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GoMp) G^,iib) + G^,12(P) + Gt,Ap) (GoMp) + GoMp)) , (B5) 



d^GoMp) 



dQ 
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and 
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2m 
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G'o,22(p) G^,22b) + Gi^^M) + G'^,12(P) (G'o,ii(p) + Go,22b)) . (B7) 



With these identities and substituting the BCS gap equation given by Eq. into Eq. (|^. 
we find 
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G'o,ii(A; + g)G'o,ii(A: + q) + Go,i2(A; + g)Go,i2(A; + g) Go,22(fc), (B8) 
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G'o,i2(fc + g)G'o,ii(A; + g) + G'o,i2(A; + g)Go,22(A; + g) Go,i2(A;), (B9) 
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2E ^|g'o,ii(A:)Go,i2(A:) [Go,ii(A;) + Go,22(A:)] 



+2E 2 
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In deriving Eqs. ()B8|) - ()B11|) . we have made use of the following identities for BCS Green's 
functions: 



and 



Go.iil-A:) = -G'o,22(fc) 



G^,i2{-k) = GoMk)- 



(B12) 



(B13) 



Taken together, Eqs. dHU and (TBSll- dnTIl) give an explicit expression for the normal fluid 
density due to Bose fluctuations in terms of products of BCS Green's functions. 

It is of interest to relate our result for to that obtained recently by Andrenacci et 
al. jl7|, who used a direct diagrammatic evaluation of the response function definition of 
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the normal fluid density (see Appendix A). Eqs. (14) and (25) in Ref. |17| give the following 
expression for the fluctuation contribution to the normal fluid density: 



pIal = -mxtz^Q = 0) 



(B14) 



where 



xt'AQ = 0) 



1 1 



(2mVB3T. E E (2^^ + 2g.)(2A;; + 2g,)r,v(g)r(g)."j"(g) 

xGo,ii'{k + q)GQ^iiii{k + q)Go,i"i'{-k)Goj>j{k' + q)Gojj"{k' + q) 
xGo,j'j"{-k'). (B15) 



This result is based on the Aslamazov-Larkin-type (AL) diagrammatic contributions to the 
transverse current correlation function Xz,z- Using our notation as defined in the text of this 
paper, the vertex functions in Eq. ()B15jl are 

M22iq) 
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and 



ri2(g)=r2i(g) = /5 



Muiq) 



(B16) 



(B17) 



detM 

To facilitate comparison with our results for p^, we expand Eqs. ()B14|) and ()B15|) to give 
1 
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+ 4ri2(g)r22(g) [GoMk + q)GoMk + q) + GoM^ + q)GoMk + q)\ GoA-k) 

X Go,ii(A;' + g)Go,i2(A;' + q) + GqA^' + q)GoMk' + q)]GoM-k') 
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GoMk' + g)Go,22(A;' + q) + GoM^' + q)Go,i2{k' + q) 



G, 



0,22 1 



-k') 



(B18) 



Comparing with our results in Eqs. (jBljl and ()B8jl - ()BlHl . it is apparent that the expression 
for PnAL Eq. ()B18|) does not contain terms analogous to those in Eq. ()B1|) arising from 
second-order derivatives of Mjj with respect to Q. Only the terms in Eq. ()B1|) that involve 
products of first-order derivatives of Mij correspond to the AL terms in Eq. ()B18|) . Separat- 
ing the contributions from first- and second-order derivatives with respect to Q in Eq. (jBljl . 
we use Eqs. fjB16jl and (jB17jl to re-write Eq. (jBljl as 



Pn 



2m , , 



rii(g)rn(g) 



11 



dMu dM: 
dQ dQ 



4rn(g)ri2(g) 



12 



+2rn(g)r22(g)^^ + 2ri2(g)ri2(g) 

dQ dQ 



dQ dQ 

dMu dM22 
dQ dQ 



dM,2dM,2 . . .dM22dM,2 ^ 
+2Ti2[q)Tl2[q)^r^ 4ri2(g)r22(g)^:T7^ 777^ + 



-hr22(g)r22(g) 



dQ dQ 

dM22 dM22 

dQ dQ 



2m2 



E 



9 L 



dQ dQ 
1 



X I M22 + Mii- 



2M 



92 M 



(MnM22 - M12M12) 




V dQ^ ' ^'^^^ -^12 5g2 yj_ • (B19) 

Using Eqs. (jB8jl and (jB9jl . as well as Eqs. (jB12jl and (jB13jl . and taking Q = z, one can show 
that the terms in Eq. (jBlQjl that involve a product of first-order Q-derivatives of the matrix 
elements Mij are precisely equal to twice the AL expression for p^^^^ given by Eq. (IBl^ 



To summarize, we have shown that the expression in Eq. (jBH) reduces to Eq. (jBTojl . This 
is equivalent to 

1 



^Pn,AL p 



X M22 



92 M 



11 



(M11M22 - M12M12) 



12 



(B20) 



where at is given by Eq. (jB18jl . 

In Ref. [ill , it is shown that the normal fluid density is indeed twice the AL contribution 
given by Eq. (jB14jl . owing to an additional AL-type diagram that is topologically non- 
equivalent to the diagram that gives rise to Eq. (jBlSjl (see Fig. 2 in Ref. |l7|). This extra 
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diagram gives rise to a contribution to the transverse current correlation function that is 
equal to Eq. ()B15I) for the case of a contact interaction potential. Hence, the final result 
obtained in Ref. [itI] for the Bose contribution to the transverse current correlation function 
in the BCS-BEC crossover is given by = 2p^^j^. 

However, we see from Eq. ()B20|) that in addition to the AL-type contribution, our expres- 
sion for the normal fluid density includes terms which arise from second-order derivatives 
of the matrix elements of the inverse Gaussian fluctuation propagator with respect to the 
superfluid velocity Q = MtK. In the BEG limit, one finds that these contributions vanish, 
which explains why Ref. [Hf also obtains the Landau formula in Eq. (|7^ in the BEG limit. 
However, the extra terms in Eq. ()B20|) are important in the unitarity and BGS regions. Our 
numerical results for p„ which are discussed in Ref. include the contributions from all 
terms in Eq. ()B20|) . These extra terms are given explicitly by Eqs. flBTn|l and (IbTII) . It 
would be interesting to understand which diagrams give rise to these extra contributions. 
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